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Abstract. We report observations of the expected main S-bearing species (SO, SO2 and H2S) in the low-mass 
star forming region L1689N. We obtained large scale (~ 300"x200") maps of several transitions from these 
molecules with the goal to study the sulphur chemistry, i.e. how the relative abundances change in the different 
physical conditions found in L1689N. We identified eight interesting regions, where we carried out a quantitative 
comparative study: the molecular cloud (as reference position), five shocked regions caused by the interaction 
of the molecular outflows with the cloud, and the two protostars IRAS 16293-2422 and 16293E. In the cloud we 
carefully computed the gas temperature and density by means of a non-LTE LVG code, while in other regions we 
used previous results. We hence derived the column density of SO, SO2 and II2S, together with SiO and H2CO 
- which were observed previously - and their relevant abundance ratios. We find that SiO is the molecule that 
shows the largest abundance variations in the shocked regions, whereas S-bearing molecules show more moderate 
variations. Remarkably, the region of the brightest SiO emission in L1689N is undetected in SO2, H2S and II2CO 
and only marginally detected in SO. In the other weaker SiO shocks, SO2 is enhanced with respect to SO. We 
propose a schema in which the different molecular ratios correspond to different ages of the shocks. Finally, we 
find that SO, SO2 and II2S have significant abundance jumps in the inner hot core of IRAS16293-2422 and discuss 
the implications of the measured abundances. 
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1. Introduction 

Low mass star forming regions are composed by at 
' least three main ingredients: the molecular cloud from 
which protostars are born, the protostars themselves, 
and the shocked regions at the interface between the 
cloud and the outflows emanating from the protostars. 
These three regions have very different physical condi- 
tions, where temperature, d ensity and also chemica l abun- 
dances greatly differ (e.g. Ivan Dishoeck fc Blakj|l99^ . 
This paper focuses on the abundance changes occur- 
ring to the S-bearing molecules and the relevant sul- 
phur chemistry. Depending on the physical condition 
of the gas, it is believed that different types of reac- 
tions play a role in the formation of sulphur-bearing 
molecules. In molecular clouds, ion- molecule reactions 
are the most important llOppenheimer fc palgarndll974l 
IPrasad fc Huntresdll9'82HMillar fc Herbstlll99(]|) . whereas 
in the warm gas of the hot cores and shocks, neutral- 
neutral reactions play the major role in forming sulphur 
species l)Pineau Pes Forets etall Il993l ICharnlevI Il997t 



iHatchell et all Il998t IXeane et al.l l200l|) . Specifically, in 
warm gas, the abundances of H2S, SO and SO2 are sup- 
posed to increase significantly. T his is the reason why they 
are o f ten used to trace sho c ks JPineau Pes Forets et al 
1993t IChernin etaD Il994t Bachil ler fc Perez Gutierrez 
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19971) . And because of the relatively fast evolution of their 
chemistry, on time scale of tens of thousand years, they 
are good candidates to be chemical clock s to study the 
evolution of outflows llBac hiller et al.l l200l|) and hot cores 
llCharnlevlll997t IHatchell et al.lll998|) . OveraU, it is widely 
accepted that in star forming regions the formation of S- 
bearing molecules is largely determined during the cold 
collapse phase, when atomic sulphur freezes out on grains 
and probably forms H2S. When the protostar starts to 
heat its environment, H2S evaporates and it reacts with 
hydrogen atoms to give sulphur atoms. S rapidly reacts 
with OH and O2 to fo rm S O, that in turn gives SO2 by 
reacting with OH (e.g. ICharnlev 199 71. 

In this paper, we present large scale maps of several 
transitions of SO, SO2 and H2S in the molecular cloud 
L1689N, a molecular cloud located in t he p Ophiuchi 
cloud complex at 120 pc from the Sun l|Knude fc Hod 
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Il998l). Based on atomic oxygen observations, ICaux et al.l 
l|1999(l found that the gas temperature in this cloud is 
(26 ± 0.5) K and the H2 density is larger than 3 xlO^ 
cm~^. L1689N harbors two young protostellar sources. 
The first one is IRAS16293-2422 (hereinafter IRAS16293), 
a Class protostar (15 L(t,) stil l in the accretion phase 
fWalker et al.' '198^": ZhoJ Il995t INaravanan et alJ Il99at 
fCeccarelli et al.i ,2J]£)0a) . Like many other young proto- 
stars, IRAS1629 3 is a binary system with a total mass 
around 1.1 Mg ()Loonev et al.ll2000|) . whose two sources 
are separated by 5", namely a projected separation of 600 
AU. The structure of the envelope surrounding IRAS16293 
has been reconstructed based on multi frequency 
SiO, O and H2CO line observations ijCeccarelli et al.l 
l2000allJ) . In the outer region (r > 150AU), the envelope 
gas shows molecular abundances typical of cold molec- 
ular clouds. In the inner region (r < 150AU, i.e. about 
2" in diameter) the abundances of II2O, SiO and H2CO 
jump to abundances typical of the hot cores around mas- 
siv e protostars. This str ucture has been recently confirmed 
bv lSchoier et alJ l)2002|) . who modeled the continuum and 
the line emission from several other molecules. The sec- 
ond protostar, 16293E, is a recently discovered low mass 
and very young Class source sit uated South-East o f 
IRAS 16293. It was detected first bv lMizuno et all l)l990fl 
as a strong NH3 peak emission. Its protostellar nature is 
discussed in lCastets et al.l l)200lL hereinafter CCLCLOl). 

The L1689N region is complex a nd has long been 
known to house mult i ple outfiows jFuku et al.l fl986t 
Wootten LoT^ri987lEmo et alJIl 99nHHirano et a,lJ 



200lV The recent work by CCLCLOl claims that the two 
protostars IRAS16293 and 16293E drive three bipolar out- 
flows. Two of them originate from each of the two com- 
ponents of IRAS16293, while the third outflow probably 
emanates from 16293E (Fig.^. In the present article we 
will adopt the scheme outlined in CCLCLOl (Fig. 1), but 
our main conclusions are substantially unaffected by the 
actual evolutionary s tage (pre-stellar or protostellar) of 
16293E, questioned in lLis et al]l)2002|) . What is important 
in the following discussion is the presence in L1689N of 
at least one protostar, IRAS16293, and six regions which 
shows an enhancement of SiO and/or II2CO emission, and 
that are sites of shocked gas marked as El, E2, HEl^, 
HE2, Wl and W2 (Fig. HI). In particular, CCLCLOl found 
that the brightest site of SiO emission, E2, does not show 
up any H2CO enhanced emission, whereas the brightest 
H2CO emission site. El, is also accompanied by strong 
SiO emission. HE2 and HEl represent a third class, for 
only H2CO emission is detected there and no SiO. The 
goal of the present work is to study how the abundances 
of S-bearing molecules change in all these sites, compared 
with the abundances in the IRAS 16293 protostar and in 
the cloud. 

The article is organized as follows. The observations 
are presented in Sect. 2, the results are presented in Sect. 
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^ This region is not considered farther in this work because 
too weak. 



Fig. 1. Sketch of the region seen face on. The dark grey 
outflows have been clearly identifled by CCLCLOl. Light 
gray used for the outflow emanating from IRAS16293 and 
the W2 source indicates that these are only assumptions. 



3, and the column density determinations are detailed in 
Sect. 4. In Sect. 5 we discuss the results, i.e. the observed 
changes of SO, SO2 and H2S abundances with respect to 
the previously measured SiO and H2CO abundances and 
what this may teach us. 

2. Observations 

Observations of large scale maps and speciflc positions in- 
side the L1689N region were performed with the IRAM 
and SEST telescopes. We measured the emission of the 
following molecules: SO (82 2i and 65 54 tran- 
sitions), SO2 (5i^5 — > 4o,4 transition) and H2S (22,0 
2i^i and li_o — ^ lo,i transitions). We also observed the 
•^^SO 32 2iline in a few positions in order to estimate 
the opacity of the main isotopic line. The coordinates 
(Aa, AS) of all maps shown here are offsets relative to 
the position of the 16293B component of the binary sys- 
tem IRAS16293 at a(2000.0) = 16'^32'"22".6, (5(2000.0) 
= -24"28'33" (the 16293A componen t is located 4" Sout h 
and 2" East from the B component - iLoonev et alJEoOfll) . 

In June 1997 we obtained a map covering 120" x 100" 
in the SO2 molecular line with the IRAM-30m telescope. 
In November 2001, May and September 2002 we per- 
formed additional IRAM observations of the SO2, '^^SO 
and H2S li,o — > lo,i lines at some specific positions in 
L1689N, namely the shocked regions El, E2, Wl, W2 and 
IIE2 (see Fig. 1 and the Introduction), the two protostars 
IRAS16293 and 16293E, and the molecular cloud (at the 
position Aa = 120", AS = 0"). In the following we will 
refer to these eight positions as the "key" positions. The 
SEST telescope was used in July 1998 to map an area cov- 
ering 300" X 200" in the SO molecular lines and to map 
a smaller area covering 120" x 100" in the H2S molecu- 
lar transition. Table Q summarizes the different observed 
molecular transitions, together with their frequencies and 
upper level energies. 

All data were obtained in the position switching mode 
with an OFF position located at Aa = -180", AS = 0" 
from the center of IRAS 16293. Using the frequency switch- 
ing mode we checked that this position is free of any C^^O 
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and ^^CO emission. Even though our observations are not 
Nyquist sampled, we smoothed our data to the largest 
beamsize when necessary, which is accurate enough, ex- 
cept in the case of pathological source morphologies. Both 
with the IRAM and SEST telescopes, pointing and focus 
was checked every 2 hours, using planets, maser sources 
and strong extragalactic continuum sources. The point- 
ing corrections were found to be always smaller than 3" 
and 5" with the IRAM and SEST telescopes respectively. 
Polynomial baselines of order 3 or less have been sub- 
tracted from the spectra. The observing parameters for 
both telescopes are listed in Tabled In the following, all 
intensities will be given in main-beam brightness temper- 
atures. Below we give some information which are specific 
to the IRAM-30m and SEST observations. 
IRAM Observations: The IRAM 30-meter telescope 
is located at an altitude of 2920 meters near the sum- 
mit of Pico Vclcta in Southern Spain. The '^"'SO 32 — > 2i, 
SO2 and H2S molecular emissions were observed simul- 
taneously with the 3, 2 and 1mm SIS receivers respec- 
tively, available at the IRAM-30m. The SO2 emission ob- 
served in June 1997 was observed simultaneously with 
other molecular transitions not presented here. The image 
sideband rejections of all receivers were always higher than 
10 dB. Typical system temperatures were in the range 
Tsys ^ 300 — 600 K. The three receivers were connected 
to units of the autocorrelator, set to provide spectral res- 
olutions of 40, 40, and 80 kHz at 98, 136, and 217 GHz 
respectively. At those frequencies, the velocity resolutions 
are of the order of 0.1 - 0.2 km s^^. 

SEST Observations: The SEST telescope is a 15-meter 
single dish millimeter telescope operated jointly by ESQ 
and a consortium of Swedish institutions. It is located at 
an altitude of about 2450 meters in La Silla (Chile). To 
obtain the SO and H2S maps shown here, we used the dual 
3/1. 3mm receiver. System temperatures were in the range 
Tsys ^ 300 — 500 K. This receiver was always connected to 
the high resolution acousto-optical spectrometer available 
at SEST, which provides 2048 channel of 43 kHz width 
giving a 86 kHz spectral resolution. At the frequencies 
considered here, that corresponds to a velocity resolution 
of 0.1 - 0.2 km s~^, comparable to that of the IRAM data. 

3. Results 

Figure El shows the velocity-integrated intensity maps of 
the SO 32 ^ 2i, SO 65 ^ 54 and SO2 5i,5 ^ 4o,4 lines. 
The spectra (not smoothed) of the six transitions of 
TableQlobserved towards the "key" positions are displayed 
in Fig. 131 and the relevant line parameters are reported 
in Table [21 The SiO 2^1 observations, published by 
CCLCLOl, have been added on Fig.Oand in TableElfor 
comparisons. 

The main-beam brightness peak temperature Tmb, 
and the Hnewidth Av (FWHM) of Table 12 were estimated 
from a gaussian fit to the whole line profile, even in the 
case where profiles are double-peaked. We also estimated 



the width of the equivalent gaussian with the same inte- 
grated and peak intensity as the observed line and found 
that the two methods differ by less than 30% with the ex- 
ception of the SO 65 — > 54 line in Wl which has an equiv- 
alent width of 2.5 km s^^ (versus 1.3 km s~^ in Tablej^l). 
On the contrary, the velocity-integrated intensity is the 
area integrated in the interval of velocities at the "zero- 
intensity" level. The velocity vlsr is the velocity at the 
peak or at the absorption dip for double-peaked profiles. 
They are roughly the same for each observed transition 
and position, with the exception of E2(HVC), whose emis- 
sion is red-shifted by ~ 6 km s^^. 

The first remarkable result is that, against probably 
naive expectations, the shocked regions are not necessar- 
ily associated with evident enhancements of the emission 
of the observed SO, SO2 and H2S lines. For example, 
the E2 region, the brightest SiO 2-1 emission peak (see 
CCLCLOl) is undetected in SO2, and only marginally de- 
tected in the lowest SO transition. 

The second result is that IRAS16293 is the only site 
where all the molecular lines have been detected, including 
the high-energy transition H2S 22,0 '^i,! - Quite interest- 
ingly, we find that not only the line profiles are broader 
than in the cloud but their widths tend to increase with 
the upper energy level of the transition. We will discuss 
in the §5 the implication of this behavior. 

We now discuss in more detail the results toward the 
eight "key" positions. 

3.1. SO emission 

The SO 32 2i emission is relatively bright (~ 3K) all 
over the cloud (Figs. |21and|2Il. All the line profiles are 
relatively narrow (with linewidths ranging between 1 and 
2kms^^) and peak around the systemic cloud velocity 
at Vlsr — 3.9kms~^. The absorption dip around the sys- 
temic velocity seen towards IRAS 16293 is probably due 
to the absorption by f oreground material, as also seen in 
the H2CO transitions ijLoinard et al.ll200l|) . On the other 
hand, the absorption dip seen in the ^''SO towards Wl 
may point to a very large SO column density in that di- 
rection or two different kinematic components. The con- 
trast of the velocity-integrated SO molecular emission at 
the various positions in the cloud is relatively low: the 
emission peaks towards IRAS16293 but otherwise it is 
widespread in the region between Wl and E2 (Fig. 0). 
Remarkably, there is not a much brighter SO 32 — > 2i 
emission towards the shocked regions with respect to the 
molecular cloud. The only two clear exceptions are rep- 
resented by El and E2, where a high velocity component 
appears in the low energy SO transition (Fig.jJJ. In the fol- 
lowing we will analyze separately the E2 low velocity com- 
ponent centered at 3.9kms~^ (hereinafter E2(LVC)), and 
the E2 high velocity component at lOkms"^ (hereinafter 
E2(HVC)), because the two components can be clearly 
disentangled. On the contrary, in El the situation is less 
clear and hence we didn't pursue a separated analysis of 
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Fig. 2. Integrated intensity maps of SO 82 ^ 2i(a), SO 65 — > 54(b) and SO2 61,5 — > 4o,4(c). (a) and (b) : first level 
is 0.5 K km with level step of 1.0 K km s~^. (c) : first level is 0.3 K km s~^ with level step of 1.0 K km s~^. 
In each map, the black points represent the observed positions. The two stars symbols show the position of the two 
protostellar sources. The arrows point to the eight "key" positions, whose spectra are displayed in Fig. |3| 
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Transitions 


SO 32 2i 


SO 65 54 


^*S0 32 -> 2i 


SO2 5l,5 


4o,4 


H2S 22,0 


-^2i,i 


H2S li,o — > lo,i 


V (GHz) 


99.299 


219.949 


97.715 


135.696 




216.710 




168.762 


E„j,/k (K) 


9 


35 


9 


15.6 




84 




24 


ncr (cm"^) 


10^ 


10*^ 


lO'^ 


2 X 10*^ 




10® 




2 x 10*^ 


Telescope 


SEST 


SEST 


IRAM 


IRAM 


IRAM 


SEST 


IRAM 


IRAM 


T.,. (K) 


160 


250 


140 


400 


270 


220 


600 


700 


Beam (") 


51 


24 


26 


18 


18 


25 


11 


14.5 


Befl 


0.75 


0.5 


0.75 


0.59 


0.71 


0.5 


0.57 


0.65 


Ai/ (kHz) 


86 


86 


40 


39 


39 


86 


42 


42 


Aq X A<5 (" X ") 


280 X 190 


288 X 216 




120x100 




120x100 






sampling (") 


24 


24 




12 




24 







Table 1. Observed molecules together with the observing parameters. The critical densities (ncr) are computed for 
temperatures between 50 and 100 K using the coUisional coefficients quoted in Sect. 4. The area mapped (Aa x Ai5) 
in each transition and the sampling of the map are reported in the last two lines. Boff is the beam efficiency of the 
telescope. The symbol "-" implies that only observations towards the "key" positions have been obtained. 



the two components. We emphasize that the E2(HVC) 
component is detected at more than 3 RMS all along the 
red lobe of the outflow and it is also observed in the SiO 
transitions (CCLCLOl). This high velocity component is 
not seen in SO 65 84 probably because of unfavorable 
excitation conditions. Indeed, we will show in the next sec- 
tion that the gas density in E2(HVC) is lower than in the 
ambient gas (see Table This fact, already recognized 
by CCLCLOl, leads to the supposition that the E2(HVC) 
may represent the wind shock, whereas the lower velocity 
component, E2(LV C), may trace the so called cloud shock 
l)Hollenbachlll998l) . Finally, we also see hints of red and 
blue wings originating from the molecular outflows driven 
either by IRAS16293 or 16293E in the direction of Wl, 
El and HE2 whereas no wing is seen in W2. 

In the higher energy line SO 65 — s- 84 the emission is 
peaked towards IRAS16293 and Wl (cf. Fig.EJ while the 
emission from the cloud is, as expected, weak. Relatively 
strong emission is also detected around El. Moreover, the 
SO 65 — > 54 spectra in those positions have pronounced 
red and blue wings. The same happens in HE2 where the 
SO 65 54 spectrum shows clearly the presence of a blue 
wing. The emitting region around IRAS16293 has a char- 
acteristic size of a X 5 = 39" x 16" namely it is not resolved 
in (5. Finally, the observed SO spectra towards 16293E 
are similar to those in the cloud: narrow, intense in the 
SO 32 — > 2i transition and weak in the SO 65 84 tran- 
sition, and no signs of wings are evident. Overall, 16293E 
does not seem to have any enhanced SO emission com- 
pared with the cloud. 

3.2. SO2 and H2S emission 

Unlike the SO 32 — > 2i , the SO2 5i,5 —> 4o,4 emission 
is limited to a few spots (see Figs. |2 and specif- 
ically toward IRAS16293, El and W2, probably be- 
cause of the excitation conditions (the critical density 
of the SO2 5i,5 — > 4o,4 transition is about 20 times the 
SO 32 2i critical density). Unfortunately, because of 
the lack of other observations around W2 we have no idea 
of its extent. Remarkably, no significant SO2 emission is 



detected towards E2. SO2 has also been detected at the 
cloud position. The linewidth (1.4 km s~^) is equivalent to 
that of N2H+ 1 ^ - a good tracer of extended, cold and 
quiescent gas - observed at the same position (CCLCLOl). 

Probably because of similar critical densities, the 
SO2 5i_5 4o,4 and SO 65 84 emission shows up in the 
same regions. In addition, with the exception of El, the 
SO2 5i.5 — > 4o.4 spectra show characteristics similar to the 
SO 65 — > 54 spectra, like for example the linewidth, in 
all "key" positions. On the contrary, in El the SO line 
profiles are characterized by a strong ambient component 
with blue- and redshifted wings of much lower brightness, 
whereas the high- velocity component appears as bright as 
the ambient component in the SO2 line, so that the whole 
SO2 5i_5 4o.4 linewidth is approximately twice that of 
the so' 65 -> 54 . 

IRAS16293 is the only position where H2S 22,0 2i_i 
has been detected. The line is relatively bright and rather 
broad {Av = 5.0kms~^). There is no emission detected 
in shocks nor in the ambient cloud. Conversely, the lower 
H2S transition, observed only toward IRAS16293, 16293E 
and the cloud, has been detected in all these positions 
with a double-peaked profile towards the protostars and 
possibly the cloud position. The absorption dip is likely 
due to the optical depth of the line. In fact, iMinh et alJ 
C199I) have observed the same H2S transition towards sev- 
eral star forming regions and found an opacity of ~10 for 
similar column densities. 



3.3. Comparison with SiO 

Contrary to the sulphur bearing species, the emission of 
SiO 2 ^ 1 is stronger in the shocked regions than to- 
wards IRAS16293. The spectra in El and Wl show high 
velocity wings much more marked than for the SO and 
SO2 lines. The high velocity component of E2 seen in 
SO 32 2i is very strong in SiO 2^1 and have sim- 
ilar Av (~ 6.5 knis~^) and vlsrI'^-" 10 kms~^) for these 
two transitions. 
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:!2 - 2i 


().-, - 54 


■.'>2 - 2i 


")i,r, — i(),4 
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2 — 1 


IRAS 
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3.7±0.1 
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16293-2422 
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Tmb 
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- 
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- 
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Av 
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- 
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- 


- 


< 0.1 
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- 
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El 


Av 


1.8±0.2 
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- 


- 
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VLSR 


4.0 


3.8 


3.9 


4.0 


- 


- 


3.6 




/ TmbSv 


7.3±1.2 


5.6±1.0 


0.8±0.2 


2.2±0.4 


< 0.3 


- 


3.3±0.3 


E2 


Tmb 


2.4±0.2 


1.1±0.2 


0.9±0.04 


<0.2 


< 0.4 


_ 


0.53±0.06 


LVC 


Av 


1.3±0.1 


0.9±0.1 


0.5±0.02 


- 


- 


- 


3.2±0.1 




VLSR 


3.9 


3.8 


3.6 


- 


- 


- 


4.0 




/ TmbSv 


3.7±0.7 


1.0±0.2 


0.5±0.1 


< 0.1 


< 0.4 


- 


1.8±0.1 


E2 




0.5±0.15 


<0.6 


<0.1 


<0.2 


<0.4 


- 


0.83±0.06 


HVC 


Av 


6.6±0.6 


- 


- 


- 


- 


- 


6.3±0.2 




VLSR 


10.4 


- 


- 


- 


- 
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/ TmbSv 


3.6±0.7 


< 0.5 


< 0.2 


< 0.3 


< 0.8 


- 


5.6±0.2 




Tmb 


2.6±0.1 


3.3±0.2 


0.2±0.04 


0.3±0.1 


<0.2 


_ 


0.3±0.1 


Wl 


Av 


1.8±0.2 


1.3±0.03 


2.2±0.2 


1.9±0.3 


- 


- 


3.3±0.8 




VLSR 


4.0 


3.3 


3.7 


3.3 


- 


- 


4.6 




/ TmbSv 


5.8±1.0 


9.0±1.5 


0.5±0.1 


0.7±0.2 


< 0.4 


- 


1.1±0.2 




Tmb 


2.4±0.2 


1.6±0.2 


0.30±0.04 


0.7±0.1 


<0.2 


- 


0.5±0.1 


W2 


Av 


1.4±0.1 


1.6±0.1 


l.liO.l 


1.2±0.1 


- 


- 


2.1±0.3 




VLSR 


3.8 


3.3 


3.9 


3.7 


- 


- 


3.5 




f TmbSv 


3.8±0.7 


2.7±0.5 


0.3±0.1 


0.9±0.2 


< 0.3 


- 


1.2±0.1 




Tmb 


3.7±0.1 


1.4±0.2 


0.6±0.04 


0.3±0.04 


< 0.2 


_ 


_ 


HE2 


Av 


0.9±0.1 


0.9±0.1 


0.9±0.1 


1.8±0.02 


- 


- 


- 




VLSR 


3.7 


3.5 


3.5 


3.6 










/ TmbSv 


5.3±0.8 


2.7±0.5 


0.6±0.1 


0.6±0.1 


< 0.3 








Tmb 


3.1±0.2 


1.1±0.2 


0.6±0.04 


0.2±0.04 


<0.3 


0.8±0.2 




Cloud 


Av 


1.0±0.04 


0.7±0.1 


0.5±0.1 


1.4±0.2 




1.1±0.4 






VLSR 


3.8 


3.8 


3.8 


3.9 




3.4 






/ TmbSv 


3.2±0.6 


0.8±0.2 


0.3±0.1 


0.2±0.1 


< 0.3 


1.0±0.2 





Table 2. Line parameters (Tmb (K), At; (km s~^), VLSR(km s~^), / TmbSv{K km s~^)) for the various transitions 
observed at the "key" positions: the two protostars IRAS16293 and 16293E, the shocked regions El, E2, Wl, W2 and 
HE2, and the reference position in the cloud. The SiO observations were previously pubhshed in CCLCLOl. When no 
signal is detected, we give an upper limit for the intensity equal to 3 RMS of the relevant spectrum. The parameters 
are taken from the spectra not smoothed. The symbol "-" implies that the relevant transition has not been observed 
or detected. 



4. Column densities 

In this section we estimate the column densities of the 
observed species, namely SO, SO2 and H2S, as well of 
the SiO and H2CO previously observed by CCLCLOl, in 
the "key" positions. In order to do that, we first estimate 
the density and temperature of the gas in each position. 
To derive the gas temperature and density, as well as the 
column densities of the different species, we used the theo- 
retical predictions from an LVG (Large Velocity Gradient) 
model described in detail in the next paragraph. The next 
three paragraphs describe the derivation of the gas tem- 
perature and densities as well as the column density of the 
observed species in the cloud, in the shocked regions, and 



in IRAS16293 respectively. Note that we did not carry out 
the analysis on 16293E, because the observed emission in 
that position seems to be dominated by the cloud emis- 
sion. 

4.1. LVG model description 

To derive the gas temperature and density, as well as the 
column densities, we compared the observed line intensi- 
ties with the theoretical predictions from a LVG model, 
which self-consistently accounts for the excitation condi- 
tions (density and temperature) as well as the line opac- 
ities. For the escape probabiUty /3 we used the following 
function of the line optical depth: /? = i~^^pJ~^'^1 ^ valid 
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Fig. 3. Spectra (not smoothed) of SO 82 ^ 2i SO 65 -> 64 ^''SO 82 2i SO2 5i,5 ^ 4o,4 H2S 22,0 -> 2i,i H2S 
li^o Iq.i and SiO 2 ^ 1 at "key" points (see text) in the L1689N molecular cloud. The vertical dashed line 
on each spectrum shows the position of the cloud systemic velocity {^fLSR—'i-^ km s~^). 



in the case of a homogeneous and isothermal semi-infinit e 
slab. More details can be found in lCeccarelh et all l)2002|) . 
Note that the dust emission is neglected in the present 
computations. 

The spectroscopic data of SO, SO2, H2S, SiO 
and H2CO are all taken from the JPL catalogue 
(http:/ / spec, wl.nas a. gov/ ftp /pub/ catalog/ catform.html; 
[ Pickett et al.lll99^ . The c ollisional coeffici ents are from 
lOreenl ill 9941^ for S O IPalmal ill 9X7h for SOo. iTiirner et a.lJ 
lll992|) for SiO, a.nd lOreenI lll99lh for H2CO respectively. 
The collisional coefficients of H2S are not available in 
literature, and we estimated them from the collisional 
coefficients of H2O, multiplied by a factor 5 (following 



the discussion in lTurneT<IT993) . This is a very rough 
approximation, giving rise to unfortunately very rough 
estimates of the H2S column density. Furthermore, the 
collisional coefficients used for the other molecules are 
not available in the full range of temperatures probed 
by our observations. In particular the low temperature 
regime is often missing. In this case we have extrapolated 
the coefficients at the lowest available temperature with 
a VTkin law. Finally, our code considers the first 50 
rotational levels of each molecule, unless the relevant 
collisional coefficients are available for a lower number of 
levels. 
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Density 


T 


Nso 


Nso. 


Nh.s 


Nsio 


NnaCO 




(lO^cm-3) 


(K) 


(10"cm-2) 


(IQi^cm-^) 


(10"cm-2) 


(10"cm-2) 


(10"cm-2) 


IRAS16293 


250 


100 


13000 


4100 


4000 


110 


750 


El 


1.8 


150 


15 


7.3 


< 96 


0.42 


7.0 


E2(LVC) 


1.0 


150 


3.2 


< 0.44 


< 260 


0.21 


< 0.86 


E2(HVC) 


0.4 


100 


7.8 


< 11 


< 1200 


1.0 


< 1.8 


Wl 


2.6 


150 


24 


2.4 


< 72 


0.15 


3.4 


W2 


1.8 


80 


8.7 


2.9 


< 76 


0.15 


3.6 


Cloud 


0.3 


26 


30 - 80 


2 


1.5 


0.02 


2 



Table 3. Gas temperature and density, and column densities of SO, SO2, H2S, SiO and H2CO in six "key" positions. 
The column densities are beam-averaged in the shocked regions (on 27" for SiO, 26" for H2CO, 24" for SO, 18" for 
SO2 and 25" for H2S), whereas they are corrected from the beam dilution in IRAS16293. Errors on the estimates are 
around 15% if considering only the statistical uncertainties. However, a factor two of uncertainty has to be considered, 
when the uncertainty on the density and temperature of the emitting gas is considered. 



4.2. Cloud 

To estimate the gas density and temperature, and the SO 
column density in the cloud reference position we simul- 
taneously best-fitted the ^''SO and ^^SO 82 -> 2i obser- 
vations (that give a line opacity for '^^SO t ~ 1.5), and 
the -^^SO 82 — > 2i and 65 84 observations (whose ratio, 
together with the ^^SO 82 2i intensity, gives pairs of 
gas temperature and density value s). Note that we used 
the elemental rat i o ^^S/^^S = 22.5 Jwilson & Roodlll994i 
IChin et al.lll99^ | Lucas fc Lisztl Il998|) . and we assumed 
that the emission fills up the beam, when comparing ob- 
servations obtained with different telescopes. The result 
of the modeling is shown in Fig. 0] We found that the 
ensemble of our observations constrains the SO column 
density in the cloud to be in the range 8 to 8 xlO^"' 
c m~^. Adop ting the gas temperature of 26 K, as found 
bv lCaux et al . (1999). we derive a density of 8 x lO^'cm"^. 

The column densities of the other species, namely SO2 
and H2S, have been derived assuming the same density 
and temperature. The SO2 abundance was found to be 
4 X 10^*, a value equivalent to the ones obtained in other 
clouds (see Table 0J. In the case of H2S we cross-checked 
the column density derivation with LTE computations, 
because of the rough approximation of the coUisional co- 
eflticients. For this we assumed an excitation temperature 
of 5 K , as found by .Swade. (.1989.1 and Hirahara et al. 
l)l995|) in similar molecular clouds (L184N and TMC-1) for 
several molecules whose fundamental transitions have A- 
coefficicnts similar to that for the H2S li^o — > lo,i transi- 
tion. In fact lMinh et al.l l|l989l) used an excitation temper- 
ature of 5 K for the H2S li_o — > lo,i transition to compute 
the column density of H2S in the molecular clouds L184N 
and TMC-1, which have temperatures and densities sim- 
ilar to L1689N. Note that increasing the excitation tem- 
perature by a factor two, would decrease the H2S column 
density by the same amount. We assumed a H2S ortho-to- 
para ratio equal to 8. Finally the SiO an d H2CO column 
densitie s are taken from the modelin g bv lCecca rclli et al] 
l)2000bj) and ICeccarelli et al.l l)200l|) respectively 2. The 

Alth ough no SiO emis s ion wa s detected towards the cloud 
position, ICeccarelli et al ] J2000bl) modeled the observed SiO 



Nso=3xlO"cra-2 Nso=8xlO"cm-2 



4.0 4,5 5.0 5.5 6,0 4.0 -1.5 5.0 5.5 G.O 

Log(I)enaitj [10*cm-']) 

Fig. 4. Cloud reference position modeling. Theoretical 
predictions as function of the gas temperature and den- 
sity for two different SO column densities. The curves re- 
fer the values of three observed quantities taking into ac- 
count the uncertainties. Solid lines: observed integrated 
intensity ratio SO 82 — > 2i over SO 65 ^ 54 , between 8.2 
and 5.2. Dashed lines: escape probability, derived from the 
34SO and ^"^SO 82 ^ 2i line ratio, between 0.1 and 0.8. 
Dot-dashed line: the minimum integrated intensity of the 
SO 82 2i transition, 2.6 K km s^^. Note that the in- 
tensity increases with increasing temperature, i.e. towards 
the upper part of the panel. Analogously, the lower limit of 
the escape probability (0.1) in the left panel is not shown 
because it lies outside the plot. 



(from Jup = 1 to 8) emission towards IRAS16293, taking into 
account the envelope physical structure, and found that the 
SiO abundance is 4.0x10"^^ in the outer envelope. In analogy 
with other molecules for which the abundances found in the 
outer envelope are equal to the (measured) abundances in the 
surrounding molecular clouds, we assumed 4.0x10"^^ for the 
SiO abundance in L1689N. In a similar way, although weak 
formaldehyde emission is detected at the cloud position, we 
used the abundance found in the IRAS16293 outer envelope 
for the cloud too. 
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Cloud 



L1689N 
L134N 
TMC-1 



x(SO) 
(10-^) 



X(S02) 

(10-3) 



X(H2S) 

(10-3) 



x(SiO) 
(10-12) 



x(H2CO) 

(io-») 



the column density of this molecule is very poorly con- 
strained. 



6 - 16 
0.6 - 10 
0.3 - 4 



0.4 
0.3 - 2.5 

2 - 6 



0.3 

3 

0.7 



4 

< 3.6 

< 2.4 



0.04 

2 
7 



Table 4. Abundances in molecular clouds. References: 
L1689N: SO, SO2 and H2S from this work, SiO from 
ICeccarelli etall l|2000bj) and H2CO from [Ceccarcll i ct al 



SO, 



abu n dance s from ISwade 



200J; L134N: SO and 

19891), H:, S from 'M inh et all ('1989') , SiO from 
Fziurvs et alJ lll989). H2CO from .Ohishi et all 1^ 992): 
TCM-1: SO abunda nces from iHirahara et al JJ1995) and 
IPratap et al l lll997l). SO2 from Ijrvine et a ll lll983i). H9S 
from Minh et alJ lll98 9). SiO from Ziurvs et alJ |l989^1 . 
H2CO from lOhishi fcKa,iful l|l99Sh . 



column densities derived for all species are reported in 
Table 131 and the abundances in Table 01 

4.3. Shocked regions 

We have used the gas temperature and density derived 
by CCLCLOl in El, E2(LVC), E2(HVC), Wl, W2\ i.e. 
in the shocked regions, by comparing the SiO observed 
lines (J from 2 to 5) with the relevant LVG model. The 
uncertainty in the derived densities and temperatures is 
around a factor two. We did not try to derive indepen- 
dently the gas temperature and density for each molecule, 
because of the lack of enough usable transitions. In fact, 
of the two SO transitions, the lowest lying transition is 
dominated by the cloud emission so it cannot be used to 
probe the shocked gas, and only one transition of SO2, 
H2S and H2CO respectively has been detected. We there- 
fore derived the column densities of those last four species 
assuming the same density and temperature derived from 
the SiO observations. This is an approximation that does 
not take into account the possible structure of the shocked 
gas, but the results are indeed not much affected by this 
assumption. In principle each molecular species may orig- 
inate in slightly different physical conditions and hence 
the computed column density may be consequently mis- 
evaluated. In practice, though, the error associated with 
this approximation is lower than about a factor two. For 
example. Lis et al. (2002), using H2CO transitions, found 
different values for the temperature (45 K) and density 
(2 X 10^ cm-^) in El. Even taking those values, the de- 
rived column density of the four molecules listed in Table 
El wouldn't change by more than a factor two with respect 
to those quoted in the table. Finally, we also assumed that 
all the used lines are optically thin (e.g. .Bla ke ct al. 1994). 
The column densities, averaged on the relevant beam, de- 
rived for all species in each of the four shocked sites are 
reported in Table Unfortunately, since we don't have 
observations of the H2S lowest transitions in the shocks. 



_4.4. IRAS 16293 

As mentioned in the Introduction, the envelope of 
IRAS16293 is formed by (at least) two components: an in- 
ner core characterized by a high temperature and an outer 
cold envelope with a temperature close to the cloud tem- 
perature. The hot core is a small region (~ 2" in diameter) 
whose te mperature is about 100 K and the density 2 .5 x 
lO^cm-3 llCeccareUi et alJl2000at fSchoier et al.ll2002() . At 
that temperature the grain mantles eva porate injectin; 
into the gas phase their constituents llCeccarelli et al 



^ Since no SiO emission is detected in HE2, the density and 
temperature cannot be constrained in this position. 



l2000albl ). The outer envelope is more extended (~ 20") 
and colder, with a temperature of 30 K (a little larger than 
the L1689N temperature) and a density of 2.5 x 10^ cm-'^ 
ijCeccarelh et al.ll2000at fSchoier et al.ll2002|) . The molecu- 
lar abundances in the outer envelope are similar to those 
in the molecular cloud, with the possible exception of the 
deuterated molecules. On the contrary, in the inner core, 
several molecules, those believed to be released from the 
grain mantles or formed from those evaporated molecu- 
les, undergo a jump in their abundances. In the follow- 
ing we focus on the abun dance of the S-bea ring molecu- 
les in the inner hot core. iBlake et al.l ljl994l) . using mul- 
tifrequency observations, constructed the rotational dia- 
grams of SO and SO2 and derived rotational tempera- 
tures relatively large (~80 K). It is therefore very likely 
that the bulk of the emission for both molecules origi- 
nates in the inner core. Furthermore, all lincwidths are 
larger than 3 km and increase with the upper level 
energy of the transition, once again arguing for an in- 
ner core origin. Pursuing this hypothesis we computed 
the column densities using these rotational diagrams, cor- 
rected for the beam dilution, reported in Table In es- 
timating the H2S column density we assumed that the 
22,0 ~^ 2i.i line emission originates entirely in the hot 
core, as it is suggested by its observed relatively large 
linewidt h. Using then the estim ates of the H2 column den- 
sity in lCeccarelli et alJ l)2000a|) (7.5 x lO^^cm^^) we derive 
the abundances reported in Table O Note that, prov ided 
that we use our volume density in the Schoier ct aL (20021 ) 
model, our derived abundance compare extremely well 
with their abundances (which, as a result, is multiplied 
by a fac tor five), supp orting the validity of our method. 
Actually ISchoier et al.l found abundances a factor 5 lower 
than ours, because they used a factor five larger den- 
sity than us in the inner core. The two density estimates 
differ because of the different diagnostics used to derive 
them: ISchoier et alJ used the continu um spectral energy 
distribution, whereas ICeccarelli et alJ used the water line 
spectrum. When CO lines are used instead. FSchoier et alJ 
found the same density in the inner region (Schoier pri- 
vate communication). W e use therefore the estimate by 
ICeccarelh et~all l|2000a|) . but keep in mind that we may 
be overestimating the abundances by a factor 5. Finally, 
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SO/H2 


SO2/H2 


H2S/H2 


SiO/H2 


H2CO/H2 


IRAS 16293 


1.7 X 10"" 


5.4 X 10" ' 


5.3 X 10" ' 


1.5 X 10"** 


1.0 X 10"' 


LI 157- mm 


5.0 X lO"'' 


3.0 X 10"* 


1.1 X 10"* 


_ 


4-6 X 10"^ 


Orion-KL 


1.5 X lO"'^ 


9.4 X 10"* 


5.0 X lO"'' 


6.0 X 10"^ 


7.0 X 10"^ 


G10.47 


3.0 X 10~^ 


1.8 X 10~* 


> 4.0 X 10"'-* 






G29.96 


> 5.0 X 10"^ 


2.0 X 10"* 


> 1.5 X 10"* 






G75.78 


3.0 X 10"^ 


> 2.0 X lO"'' 


> 1.0 X lO"'' 






G9.62 


> 4.0 X 10"" 


1.0 X 10"* 


> 8.0 X lO"'' 






G12.21 


> 8.0 X 10~^° 


> 5.0 X 10"^° 








G31.41 


4.0 X 10"^ 


1.2 X 10~* 


> 3.0 X 10"^ 






G34.26 


6.0 X 10"^ 


1.5 X 10"* 


> 2.0 X lO"'' 






BF- Class 


3.0 X 10"^ 


4.0 X 10"^° 


2.0 X 10"^ 







Table 5. Abundances of SO, SO2, H2S, SiO and H2CO with respect to H2 in the hot core of IRAS16293. The 
abundances of SO, SQ^ and H^S are estimate d by the present work observations, assuming a H-:> column den- 
sity of 7.5 X lO^^cm"^ llCeccareHi et al.ll2nnna^ . The SiO and H2CO abundances are taken from ICeccarelli et alJ 
ll20^0a^ and ICeccarelli et alJ ll2non iJ) respectively. References of the other sources : L1157-mm(I RAS 20386-1-6751): 
Bachiller fc Perez Gutierre3l|l997|) : Orion-KL: SO, SO2, SiO and H2CO from ISutton et al . (1995), H2S from 



M 



0: 



nh et al.l tmd) : G10.47. G29.96, G75.78, G9.62, G12.21, G31.41 and G34.26: 



Buckle fc Fiilierl ll2n03 



Hatchell et al.. (4998,1 : BF-Class 



the SiO and HqCO column de nsities are taken from the 
ICeccarelli et aH l|2nnnalhl. EinT^ modeling. 



5. Discussion 

The abundance ratios between the observed species are 
reported in Table 0] |31 and El They show variations up to 
two orders of magnitude. In the following we analyze in 
detail the variations associated with the cloud, the shocked 
regions and the protostar IRAS16293 respectively. 

5.1. The cloud 

Adopt ing the H2 column density derived bv ICaux et alJ 
l(l999l) by means of atomic oxygen observations in L1689N, 
5 X 10^^ cm~^, and using the column densities derived in 
the previous section (Tab. ^ we obtain the abundances 
reported in Tab.0] In the table we also report, for compar- 
ison, the abundances measured in L134N and TMC-1, two 
among the best studied molecular clouds. With the possi- 
ble exception of formaldehyde, which seems to be under- 
abundant, L1689N has abundances typical of other cold 
molecular clouds and it is therefore interesting to study 
how those abundances change in the region either because 
of the presence of shocks or the presence of a protostar 
(IRAS16293). 

5.2. Shocked regions 

In the shocked regions, with the exception of E2, it is 
practically impossible to derive the absolute abundances 
of the observed molecules, because of the difficulty to de- 
rive reliable estimates of the H2 column densities. In fact, 
usually the H2 column densities are estimated from the 
CO millimeter observations converting the measured CO 
column densities into H2 column densities. The method re- 
lies on the capacity to disentangle the contribution of the 



cloud from the shocked gas in the CO emission. However, 
this is feasible when the shocked gas emits at relatively 
large velocities, i.e. when the two contributions can be 
separat ed based on their spectral proper ties (see for ex- 
ample iBachiller fc Perez Gut ierredl 19971) . In the specific 
case of the outflows in L1689N, the projected velocity of 
the shocks is too small and the cloud and shocked gas 
cannot be disentangled, except in the E2 position, where 
a high velocity component is present. In order therefore 
to estimate the H2 column density in the shocked regions 
we would need high spatial resolution observations, be- 
cause single dish CO measurement are totally dominated 
by the cloud emission. For this practical reason, we com- 
puted only abundances ratios, as summarized in Table El 
Note that the quoted ratios have been obtained dividing 
column densities averaged on different beam sizes (18" for 
SO2 observations and ~ 25" for the others), so that the 
SO2 ratios may be off by a factor two (in the shocked 
regions) . Fig. [3 shows a graphic representation of the de- 
rived abundances ratios in the shocked regions and in the 
cloud. 

The first result to note is that, with the exception of 
the SO2/H2CO ratio, all other ratios in the shocked re- 
gions differ by more than a factor two when compared to 
the ratios in the cloud. In particular, the largest variations, 
by up a factor 100, are shown by the ratios with the SiO, 
confirming that SiO is the best molecule to trace shocks 



(Ga,selli et aLll997tlSchilke et al.ll 997HGaTav et al.l2002ll. 

A second robust result is that SO2 is enhanced with 
respect to SO in the El, Wl and W2 shocks, when com- 
pared to the cloud (by up a factor 10). This also confirms 
a rather general theoretical prediction that SO2 is over- 
produced with respect to SO in the shocked gas, because 
of endothermic reactions that lock most of the gaseous 
sulp hur in SO'^i on timescales rela tivelv small ~ 10^ yr 
fe.g. lPineau Pes Forets et al.ll993(l . Remarkably, SO2/SO 
is lowest in the strongest SiO shock of the region, E2. 
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Region 


SO/H2CO 


SO2/SO 


SO2/H2CO 


H2S/SO 
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SiO/H2CO 
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~ 0.04 


1 
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0.01 


0.01 


El 


2 


0.5 
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0.03 


0.06 


0.06 
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> 4 


< 0.1 


- 


< 80 


0.06 


> 0.5 


> 0.2 
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> 5 


<0.1 


- 


< 150 


0.1 


> 0.9 


> 0.5 


Wl 


7 


0.1 


0.7 


< 3 


6 X 10"^ 


0.06 


0.04 


W2 


2.5 


0.3 


0.8 


< 9 


0.02 


0.05 


0.04 


L1157-B1 


0.6 - 1 


0.6 - 1 


0.4 - 1 


0.8 - 1.3 


0.2 - 0.3 


0.3 


0.1 - 0.3 


L1157-B2 


1 


1 - 3 


1 - 3 


0.7 - 2 


0.1 - 0.4 


0.1 


0.1 - 0.4 


CBS 




~1 




<1 


<1 


<1 




L1448 










0.3 






Bl 










0.1 






CepA 










0.08 






NGC 2071 










0.01 - 0.05 







Table 6. Abundance ratios in the shocked regions and in the cloud, as derived from the column densities of Tab. El 
Refere nces for other outflows: L1157-B1 and L1157-B2: Bachiller & Perez Gutnx^rrez J 1997); CB3' Codella fc Bachillen 
lll999t) : L1448, Bl and CepA: iMartin-Pintado et all Il992l) : NGC 2071: iChermn fc MassonI lfl99.'£ 




Fig. 5. The abundance ratios observed in the different shocked regions and in the molecular cloud and presented 
in Table m SO/H2CO, SO2/H2CO and SiO/HaCO on the lower panel and SO2/SO, SiO/S02 and SiO/SO on the 
upper panel. The arrows represent upper and lower limits of the SO2/SO, SiO/S02, SO/H2CO and SiO/H2CO ratios 
respectively. 



CCLCLOl argued that E2 is an older shock based on 
H2CO and SiO observations compared with models pre- 
dictions and this would agree with the low SO2/SO ratio, 
as SO2 is expected to be con verted into atomic su lphur at 
late stages > 3 x 10"* yr (e.g. lHatchell et al.lll99^ . 

The third result, is the constant SO2/H2CO ratio, 
within a factor 2, in the shocked regions El, Wl and W2, 
and in the cloud. Why SO2/H2CO has the same value in 
these shocked regions and in the cloud seems a remark- 



able coincidence for which we do not have an explanation. 
On the contrary, the constant SO2/H2CO ratio in these 
shocked regions may be due to the timescales involved 
in the formation and destruction of the two molecules, 
even though the mechanisms of SO2 and H2CO formation 
in shocks are expected to be different. On the one hand, 
formaldehyde is thought to be formed onto the grain man- 
tles, released into the gas phase because of the mantle 
sputtering in the shock, and finally converted into more 
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Time (yr) 


Formation 


Destruction 


0" 


HaS, H2CO, SiHa, SiH4 




10^ 


SO, SO2, H2S, SiO 




3 X 10* 




SO2, H2CO 


10^ 




SiO 



Table 7. Proposed schema of the formation and destruc- 
tion of SO, SO2, H2S, SiO and H2CO as function of time. 
The times are indicative and are meant to represent a 
sequence. ° Injection in the gas phase of the molecules 
sputtered from the grain mantles. 



complex molecules by ga s phase reactions on timescales 
larger than ~ 3 x 10^ yr l)Charnlev et al.lll992() . SO2 fol- 
lows a totally different route: it has been proposed that 
sulphur is released into the gas phase mainly as H2S, 
following also mantle sputtering, and then H2S is trans- 
formed into SO and SO2 via neutral-ne utral r eactions on 
timescales of ~ 10^ yr (e.g. CharnlevI Il997t) . Later on 
'--^ 3 X 10'' yr, gaseous sulphur is expected to be in atomic 
sulphur. Therefore, the constant SO2/H2CO ratio in El, 
Wl and W2 may tell us that these three shocks have sim- 
ilar ages, around 10'* yr if models are correct, whereas E2 
is a much older shock, > 3 x 10"* yr, and SO2 and H2CO 
both disappear because transformed into S and complex 
0-bearing molecules respectively. 

The overall emerging picture is that the combination of 
the variations in the abundance of different molecules can 
in first instance give an estimate of the age of the shock. 
Table [3 summarizes the situation. The times noted in the 
Table {7\ are really modeled dependent and are meant to 
represent a sequence. The simultaneous presence of SO, 
SO2, H2S, H2CO, and SiO would mark relatively young 
shocks (< 10* yr), whereas the only presence of SiO emis- 
sion would point towards relatively older shocks (> 3 x 10** 

yr)- 

Finally, comparison with other molecular outflows 
can only be very approximate, but still somewhat il- 
lustrative (Tab. inil. The only molecular outflow at 
our knowledge in which SO, SO2, SiO and H2CO 
abundances have been measured is the L1157 outflow 
ijBachfller fc Perez Gutierre j Il997|) . Looking at Table El 
L1157 Bl and B2 are rather similar to E2 with respect to 
the ratios involving SiO, i.e. they are relatively enriched 
in SiO with respect to the other shocks of L1689N and the 
molecular cloud. The likely interpretation is that E2, Bl 
and B2 are all relatively strong shocks. Yet, contrary to 
what happens in E2, H2CO, SO and SO2 are detected in 
L1157 Bl and B2, which would point towards relatively 
young shocks (Table[7I). Note though that SO seems to be 
underabundant with respect to H2CO and SiO in Bl and 
B2 when compared not only to E2, but also to El, Wl and 
W2. At the same time, SO2 is overabundant. Whether this 
is because of the stronger shock in Bl and B2 that converts 
SO into SO2 more quickly, or because the shocks in LI 157 
have a different age of El, Wl and W2, is impossible to 
say at this stage. A larger statistics on a larger number 
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Fig. 6. Theoretical predictions of the evolution of SO2/SO 
and H2S/SO ratios a s a function o f time . The curves 
are adapted from the iHatchell et alJ l)l998j) model, that 
assumes a gas temperature of 100 K and a density 
of 2 x lO^cm^^. The grey squares represent the por- 
tion of the modelling curve in agreement with ob- 
servations of the hot core of IRAS16293 (this work). 
Analogously, the stars mark the predicted values of the 
ratio s in agreement with what observe d towards L1157- 
mm (iBachiller fc Perez Gutierrezlll997l) . 

of outflow system is necessary to better understand this 
point. 

5.3. IRAS16293 

Table El quotes the absolute abundances of SO, SO2, H2S, 
H2CO and SiO in the hot core of IRAS16293 and in 
the molecular cloud. As evident from the table, the sul- 
phuretted molecule abundances exhibit a strong enhance- 
ment with respect to the molecular cloud abundances. The 
abundance of SO increases by a factor 200, SO2 by a fac- 
tor 1 300 and H2S by a factor 1700 (see also ISchoier et alJ 
I2OO2I) . 

As we mentioned in the Introduction, the sulphur 
species are expected to be useful to evaluate the age of 
hot cores. We could not resist the temptation to compare 
the measured abundance ratios with published theoreti- 
cal predictions. In Fig. Owe compare the H2S/SO and 
SO2/SO ratios obtained in the IRAS16293 hot core with 
the theoretical values predicted by chemical models. The 
observed H2S/SO and SO? /S O ratios are consistent with 
the theoretical predictions by 'Hatchell et al.' ('1998'), but 
differ substantially with those by Charnlcy (1997). Both 
model s use the sulphur chemistry described in ICharnlevI 
1I1997I) and assume that the bulk of sulphur is stored in 
iced H2S onto the grain mantles. When the dust temper- 
ature exceeds 100 K, the grain mantles evaporate releas- 
ing the H2S into the gas phase, which is slowly converted 
into SO and subsequently into SO2 by neutral-neutral re- 
acti ons. The only d ifference between the two models is 
that iHatchell et al.l use an updated rate for the reaction 
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of H2S with atomic hydrogen. This difference substan- 
tially changes the evolution of the abundances of H2S, 
SO and SO2 at times l arger than about 10"^ yr. Our ob- 
servations support the iHatchell et~al] l)l998|) model, as 
shown in Fig. The agreement between observations and 
model predictions for an age of ^ 3 x 10"* yr is remark- 
able, as independent estimates o f the age of IR AS 16293 
converge towards the same value (jCeccarelli et al . 2000a). 
On the same figure we also report the values observed 
towards L1157-mm. Based on these observations and on 
the model predictions, L1157-mm has about the same age 
of IRAS16293. We caution that the error bars are indeed 
large enough to make the agreement between observed and 
predicted values just a coincidence. In this sense, this re- 
sult may just mean an encouragement to pursue this kind 
of studies on a large sample of low mass protostars, where 
a statistical trend could be established. We have indeed 
started such a study on a sample of low mass protostars, 
whose we have derived the physical structure (i.e. density 
and te mperature profiles) and formaldehyde abundance 
profiles l)Maret et alJl2002|) . 

A direct comparison between the measured abun- 
dances in the hot cores quoted in Tab. and that of 
IRAS 16293 can only be done with respect to Orion- 
KL, as in the other sources the abundances are not 
corrected for the be am dilution (JIatchell e t al... ,1998 : 
iBuckle fc Fime3l2 003'). With respect to Orion-KL, the hot 
core of IRAS16293 is enriched in SO, SO2 and H2CO, 
whereas it is deficient in H2S. The most plausible expla- 
nation is that the different abundances reflect a differ- 
ent composition of the ices (from which II2CO and II2S 
evaporate), a fact now well documented by the study of 
deuterated mole cules in high and low mass protostars (e.g. 

However, the different gas temperature 
and density may also play a role. Finally, the SO2/SO 
and H2S/SO ratios of IRAS16293 differ by about a factor 
100 and 10 respectively with respect to the relevant ratios 
measu red in the high mass protostars of the Hatc hcU ct a l. 
l)l998|) sample. Again, the difference may be attributed to 
both different evolutionary stages and different gas densi- 
ties and temperatures. 

Following the theoretical expectations that gaseous 
sulphur in the warm gas is mainly lock ed into SO , SO2 
and II2S (e.g. ^ineau Des F orets et al.l Il993: C harnlevj 
1997 ) , while silicon is mainly locked into SiO ll|Herbst et al.l 
1989|) . the (SO-hS02-hH2S)/SiO ratio gives a measure 
of the gas phase elemental sulphur over silicon abun- 
dance ratio in the IRAS16293 hot core. This value is 
about 300, to compare with the solar elemental sulphur 
over silicon abundance ratio of 0.5. It is unlikely that we 
overestimate the overall quantity of sulphur. Conversely, 
it could have been underestimated because of a possi- 
ble optical thickness of the sulfuretted molecular lines 
and/or the presence of other important sulphur-bearing 
molecular species not considered here. Thus, the mea- 
sured (SO-|-S02-l-H2S)/SiO ratio suggests that there is 
an important deficiency of silicon in the warm gas of the 
IRAS16293 hot core. The likely explanation is that sil- 



icon is mostly depleted into the non-volatile, refractory 
cores grains, whereas sulphur is depleted mostly onto the 
vo latile grain mantles . This conclusion supports the thesis 
bv lRuffle et al.l l|l999^ that sulphur and silicon follow dif- 
ferent routes of depletion, and specifically that sulphur is 
depleted at the time of mantle formation. Another pos- 
sibility to explain this deficiency of silicon is that sili- 
con is mainly in another form than SiO in the warm 
gas such as Si, Si+ or Si02. The atomic silicon is dif- 
ficult to observe, since it has its ground transition at 
129 /im, i.e. a wavelength range obscured by the atmo- 
sphere. Observations obtained with the Long Wavelength 
Spectrometer on board ISO did not detect any signal large r 
than about 5x 10^"'^'^ erg s~^ cm~^ llCeccarelli et al 
corresponding to an upper limit for the Si column density 
of 6 X 10^^ cm~^, and an abundance lower than 8 x 10~^, 
i.e. about 500 times more than the SiO abundance found 
in IRAS16293. For Si+, it is unlikely that the UV field in 
the observed regions is strong enough to ionize the silicon. 

Finally, the total abundance of SO, SO2 and II2S in 
the IRAS16293 hot core is 2.8 x 10"^, namely more than 
ten times lower than the corresponding solar abundance 
3.4 X 10^^. Even adding up the OCS. whose abundance is 
about that of SO ijSchoier et alJl2002() the overall sulphur 
abundance is still low. Unfortunately l|Schoier et al ]|2002|) 
could not estimate the OS abundance in the hot core of 
IRAS16293, but CS is unlikely to be the main reservoir 
of sulphur. The other possible reservoir of sulphur, the 
atomic sulphur (Charnlev 1997), is extremely difficult to 
observe, because the ground state transition is at 25 /im, 
i.e. in a wavelength range obscured by the atmosphere. 

6. Conclusion 

We have presented a quantitative observational study of 
the most important S-bearing molecules, namely SO, SO2 
and H2S, in the region of L1689N. We derived the column 
density of these molecules plus SiO and II2CO molecules 
in six regions of L1689N: the cloud, the young protostar 
IRAS16293, and four shocked regions. 

We found that SiO is the molecule that shows the 
largest abundance variations in the shocked regions, 
whereas S-bearing molecules show more moderate vari- 
ations. Remarkably, the region of the brightest SiO emis- 
sion in L1689N, namely E2, is undetected in SO2, H2S 
and II2CO and only marginally detected in SO. We ar- 
gued that this is possibly due to the relatively old age 
(> 3 X 10"* yr) of this shock. 

In the other weaker SiO shocks, SO2 is enhanced with 
respect to SO, in agreement with theoretical expectations 
that predict the conversion of the gaseous sulphur mostly 
into SO2 on timescales of ~ 10"^ yr. In the same regions, 
the SO2/H2CO ratio is of order of unity. We argued that 
this may point to relatively young shocks (~ 10"* yr), 
where SO2 has already formed and H2CO has not yet 
destroyed. 

Putting together the observed combinations of the SO, 
SO2, H2CO and SiO ratios, we proposed a schema in which 
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the different molecular ratios correspond to different ages 
of the shocks. 

Finally, we found that SO, SO2 and H2S have signifi- 
cant abundance jumps (200, 1300 and 1700 respectively) 
in the inner hot core of IRAS16293. We compared the mea- 
sured abundances with theoretical models and discussed 
the derived protostar age. However, we cautioned that a 
more detailed study is necessary to draw reliable conclu- 
sions. The hot core of IRAS16293 seems to be enriched 
in SO, SO2 and H2CO with respect to Orion-KL, proba- 
bly because of a different initial composition of the ices in 
the two sources. Comparing the S0-|-S02-f H2S/SiO ra- 
tio in the hot core of IRAS16293, we found that silicon 
is largely deficient in the warm gas (by a factor ~ 600), 
supporting the thesis that silicon is depleted into the grain 
refractory cores whereas sulphur is depleted into the grain 
volatile mantles. Nonetheless, sulphur in the IRAS16293 
warm gas is also deficient. 
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